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Abstract 

We study how inhomogeneities modify large scale parameters in General Relativity. 
For a particular model, we obtain exact results: we compare an infinite string of ex- 
tremal black holes to a corresponding smooth line with the same mass and charge in 
five dimensions. We find that the effective energy density does not differ significantly. 



1 Introduction 



Less than 5% of all matter and energy in the universe is visible-the origin of the rest is a mystery. 
Dark Energy (DE) constitutes perhaps 70% of this missing mass (see, e.g., this is the energy 
needed to accelerate the whole universe outwards at its measured rate. The remaining 25 ~ 30% 
is Dark Matter (DM) (see, e.g., Q), and is needed to account for the proper motions of celestial 
bodies. 

This paper investigates how averaging out inhomogeneities in matter affects the measurement 
of mass and energy parameters. When we calculate interactions between celestial bodies, we 
generally use a smooth mass density profile in galaxies or clusters; that is, we average mass dis- 
tributions to simplify our calculations. However, Einstein's theory is nonlinear, and this averaging 
could potentially alter what energies we expect to see. The averaged stress tensor of the real metric 
need not be the same as the real stress tensor of the averaged metric; this discrepancy could be 
pertinent to the problem of DM. 

Effects of inhomogeneities and DM seem to have some characteristics in common. Smoothing 
of mass fluctuations should become less important at greater distances from the fluctuations; thus 
the effects on energy density would be greater closer to the masses, which is where halos are iden- 
tified. It also make sense that the existence of halos is contingent upon some intrinsic properties of 
mass distributions in general, because the halos are found on such different scales (galaxies, galaxy 
clusters, even superclusters). This marks a difference from DE, which is seen as a constant energy 
density throughout the whole universe, a "cosmological constant" as Einstein originally called it. 

This is not the first proposal with respect to the effects of smoothing on energy and the con- 
sequent significance for cosmology. Kolb et al. suggested that inhomogeneities could account 
for a cosmological constant and solve the problem of DE rather than DM LiJ. This proposal was 
refuted by E. Flanagan H and R. Wald and A. Ishibashi 0, on the grounds that the effects of 
such fluctuations are far too small to compensate for the huge amounts of energy missing. In the 
context of DM, F. Cooperstock and S. Tieu claimed to construct a model in which a galaxy's need 
for a halo was canceled by including some of its possible rotational energy [6|; M. Korzynski and 
others have shown this claim is erroneous \JJ. Other explorations of these issues include 

All these calculations, however, are perturbative-based upon approximations; until now, there 
has not been an exact calculation of the impact of smoothing on energy. This paper determines the 
exact energy for two scenarios, one of which is a 'smoothed' version of the other. Both scenarios 
use BPS black hole solutions because it is easy to find a solution for any sum of these special 
black holes, as the electromagnetic repulsions and gravitational attractions between them cancel 
ll9l . Specifically, we compare an infinite string of discrete charged black holes to a smooth line of 
charge, with both distributions normalized to have the same average charge per unit length. We 
work in five dimensions because such distributions are divergent in four dimensions. 
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2 The Metrics 



The solution of the coupled Einstein-Maxwell field equations which describes a set of n = I . . .N 
extremal black holes with charges qn in five dimensions^ is ifTOl : 

/3 ^ a 

as we are in five dimensions, dx^ = dx^ + dy^ + dz^ + dv^. We note that is time-independent 
and harmonic (it satisfies V^Q. = outside the source), and that, in these coordinates, the horizon 
of each black hole is shrunk down to a point. 

These are BPS solutions, with the property that gravitational attraction and electromagnetic 
repulsion between any set of them cancel [9|. This makes them time-independent: with no forces 
between them, they are in equilibrium. This is also what ensures that solutions for a set of them 
are easily found, namely by the sum in (12.11) . 

We consider two configurations: a discrete (D) and a smooth (5) solution. In (D), we have 
a discrete set of black holes positioned along a line at equal distances L from one another other, 
whereas in (5) we have smooth charged string along an axis. The continuous set can be viewed 
as a limiting case of the discrete set, namely the limit as the number of black holes per unit length 
tends to infinity and as the charge of each black hole tends towards zero, in such a way that the 
total charge per unit length remains the same. 



2.1 The Discrete Solution (D) 

For the discrete case, we choose black holes with the same charge qo', thus Q.d is given by: 

oo 

where black holes are one unit length apart on the v-axis and r is the distance to the v-axis.. We can 
restore the separation L by dimensional analysis (keeping the total charge per unit length constant): 

= ^ + r^ + {v-nLY = ^ + „L (r/L)2 + (n-[v/L])2 " ^^.3) 

Thus we can start with L = 1 and find the dependence on L simply by making the substitution 
(r,v,^Z)) {r/L,v/L,qD/L). 

To evaluate the sum in (12. 21) . we use contour integration. We find 

^ 1 /;rcot(;r(v + z>)) z;rcot(;r(v-z>)) 

1 \^-^) 



r2 + (v-n)2 2r 2r 



'in four dimensions, the solution 1.9 lOJ has the form ds^ — Q. ^dt^ + Q?dx^, Aq — Q. ' . 
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which simplifies to 

^1 n sinh(2;rr) 



r^H- (v — r cosh(2;rr) — cos(2;rv) 
Thus, using (12 .21) and (I2.5I) . the metric for the discrete solution is given by: 



(2.5) 



o 1 , sinh(2;rr) 

r cosh(2;rr) — cos(2;rv) 

2.2 The Smooth Solution {S) 

The metric for (5) is the same as the limit of the metric of {D) as the masses and the distances 
between them approach zero. From (12.31) . we find: 

as = i + \im £ , ^^^ = i+r dx , ,2 = 1 + ^- (2-7) 



As expected, the discrete harmonic function Q.d (12.61) approaches the smooth function Q5 (12.71) 
when the distance r to the v-axis is large. 

3 Energy 

Since we have normalized the two solutions so that they have the same average charge density 
along the v axis, and since the solutions are BPS, they must have the same total mass density as 
well; this does not immediately rule out a possible effect, as we discuss below. 

We check this using Weinberg's definition of total energy ifTTTl : we rewrite the metric as a 
flat metric plus a correction term h: g^y = rj^y +^/xv- We split the Einstein tensor into terms 

linear and non-linear in h: G^fy (h) + G^^^y^ (h) = T^y. Then the definition of the total (matter plus 
gravitational) energy density of a system is given by ifTDl : 

Total Energy Density = g'^^ (h) = Tqq - g'^^^ (h) . (3.1) 

We integrate Gqq^(/z) over all three directions (x,y,z) orthogonal to the line of black holes. Inte- 
grating along the entire v-axis would make the integral divergent-but this is not a problem: since 
the smooth configuration (5) is independent of v, and the discrete configuration (D) is periodically 
dependent upon it, we need only integrate over one period of configuration (D) (i.e., from to 1 
for L = 1). To perform the integral, we transform to 4-d cylindrical coordinates (r, 0, (p, and v), 
and integrate over v, 0, (p, factoring in the Jacobian r^sm{(p) for the coordinate change. Then the 
total energy per unit length along the v-axis (describing an infinite slice of space orthogonal to the 
v-axis, with a width of one unit length) can be written as the following surface integral: 

E = J dv J ded(pr^sm{(p)n\dihjj-djhij) , (3.2) 
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where H is a unit vector normal to the v-axis. Since at large values of r the integrands are the same, 
this integral yields the same result for our two cases-in agreement with the requirements of the 
BPS condition (which ensures that interaction energies cancel) and the equality of charge per unit 
length (which ensures that the energies which are not due to interactions are also equal). 

The equality of the total masses does not in itself preclude the possibility of an effect on DM: 
the energy which enters Einstein's equations, and so is relevant to DM, is the covariant matter 
energy Tqq. This quantity is not necessarily the same for our two different distributions of black 
holes; comparing them shows us how matter energy is affected by smoothing. 

3.1 Stress Tensor and Integration 

The matter energy of a system is given by the integral of Tqq over all space. In a general spacetime 
with a time-like Killing vector the following total derivative vanishes: 

= J dhV^Dp{gP^T^y^^) = J dhdp{V^gP^^T^y^^) ; (3.3) 

consequently, the contribution on a spacelike slice defines a conserved energy. In our case, since 
both the (D) and (5) metrics are static and diagonal, the total matter energy reduces to: 

£ = -y'j4^v^/Oroo . (3.4) 

The stress-tensor is usually defined with a factor j^, but as we only study ratios of energies, we 
are not interested in the overall normalization. 
We calculate Tqq using 

7)xv = F^aFypg""^ - \g^v{Fa^FySg''h^^) , F^y = d^Ay - dyA^ , (3.5) 

where A is the gauge potential (12.11) . Once again, because the metrics are static and diagonal, and 
furthermore, because the gauge potential A is purely timelike and static, this simplifies, and we 
find 

Too = FojFojg'J - ^goo{2FojFojg'^''gjj) . (3.6) 

Recalling Ao = ^Q.^^, goo = gjj = and the determinant g = —Q? from equation (12.11) . 

we find: 

E = - ld'x ^gg^ Too = 6 1 d'x . (3.7) 

Here {djSl)^ = {d^Q)^ + {dya)^ + {d,a)^ + (d^a)^. 
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3.2 Exact Energies 

As discussed above equation (I3.2I) . we integrate stress-energy over all three directions {x,y,z) or- 
thogonal to the line of black holes and integrate over one period of configuration (D) over the 
v-axis. The smooth configuration (5) is the simpler of the two; we calculate it first. As discussed 
in section 2.1, we use L = 1, and reintroduce L later by rescaling. We ignore integration over G 
and (p in both configurations because it only produces equal overall constants; we also drop the 
sin((p) factor from the Jacobian. Defining k = nq and Q5 = 1 + ^, we find (djO.s)'^ — k^/r^, and 
then (13.71 ) gives 

{r + k}'^ 

For discrete configuration (D), the calculation is more complicated. We define s = sinh(2;rr), 
c = cosh(2;rr), t = sin(2;rv), and d = cos(2;rv) and compute: 



,2_k^ [s{d - c) + 27tr{\ - cd)]^ + {iTtrst)^ 
7 {c-dY 



{djilo) = -T ■ ^^■^> 



Then, with Q.d = \ + 72^' ^'^^ measure factor r^, the energy for (D) is given by: 

Ed= dv drk^^-^ / — -y^. , . (3.10) 



Jo {c-dY{r{c-d)+ksY 

This is not easily evaluated as an indefinite integral, but we can compare Eg and Eq numerically. 
To understand the results, let us clarify what we should test. For L = 1 and various k (charge), we 
want to find the differences in energy between configurations (D) and (5) . However, a smaller k 
will certainly yield a smaller energy and smaller energy discrepancy, so to compare various ^'s we 
divide by Es because k Eg. 

Thus we study the ratio ^ /J° /J [Eq — Es)dv dr^ / ^ /q Es dv dr^ as a function of k. Different 

k represent different densities, as we have scaled the separation Lto L= 1. The results, computed 
with Maple, are consistent with zero; the actual numerical integrations are somewhat delicate 
because of divergences in the integrand (Ed diverges when r = and v = 0, 1). 



4 Physical Application 

Though we did not find an effect, we discuss how any effect that we might have found could have 
been interpreted physically. The first step would have been to choose a scale, because DM is found 
on several different levels: galaxies, galaxy clusters, superclusters. Then we would determine 
which objects are smoothed for the situation chosen. For example, to see how smoothing might 
create a halo around a galaxy, we would smooth the stars inside the galaxies, as is the common 
practice when calculating intergalactic motions. 

Our numerical results were calculated for L = 1. We can use these results, however, if we 
rescale the situation by changing the mass so that the density of the galaxy is correct. The value of 
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k in equations (13.81 ) to (13.101 ) would then be modified to provide the correct mass matching up with 
the real-life situation. The k we used was nq, and q = Gm for extremal black holes. The ratio to 
demand for similarity is fc/1 = TzGm/L. 

So, we need the mass (m) and average distance between masses (L) as input, and then we can 
see what the energy difference is for k = nGm/L. Let us apply these ideas to the situation of galaxy 
halos just discussed, testing how the smoothing of stars within the galaxies could affect the energy. 
We use planck-units, so that G = 1. The average distance between adjacent stars is about 10 light- 
years, 10^^ Planck lengths. The average mass of a star is a bit less than 10-^^ Planck masses. Our k 
is then around 10 

Next let us apply the smoothing to galaxies within galaxy clusters. The average galaxy mass 
is 10^^ Planck- masses. The average separation between galaxies is 10^ light-years, 10^^ Planck- 
lengths. So k for this situation is again 10^^^. 

A third scale at which to look is that of superclusters. Superclusters consist of about a dozen 
clusters, each with a mass of about of 10^^ Planck-masses. Usually they are arranged in strings 
with lengths of 10-100 Mpc, 10^^ or 10^^ Planck-lengths, so the distance between each is about 
10^^ or 10^^. Our is then 10"^ ^j. iq-3 

5 Conclusion 

The calculation performed in this paper is by no means flawless. Five-dimensional effects are ex- 
pected to show similar trends as their four-dimensional counterparts, but when one is discussing the 
extent of an effect, 5-d is less applicable. Also, these calculations were done with BPS solutions. 
They are easy to work with but certainly do not represent standard matter. They are maximally 
charged, for one thing, and they describe point-like masses, a physical impossibility. 

These are the first exact calculations of the effects of inhomogeneities; as such, they comple- 
ment and confirm perturbative arguments. There is room for development-our calculation was 
based on a particular model, and calculations in other models could help elucidate the exact effects 
of smoothing. 
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